Topologically Driven Swelling of a Polymer Loop 
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Numerical studies of the average size of trivially knotted polymer loops with no excluded volume 
are undertaken. Topology is identified by Alexander and Vassiliev degree 2 invariants. Probability of 
a trivial knot, average gyration radius, and probability density distributions as functions of gyration 
radius are generated for loops of up to N — 3000 segments. Gyration radii of trivially knotted 
loops are found to follow a power law similar to that of self avoiding walks consistent with earlier 
theoretical predictions. 



Although knots in polymers have been studied for sev- 
eral decades, they remain perhaps the least understood 
subject in polymer physics. Most of the work in this 
area has been directed at classification of knots, find- 
ing efficient topological invariants, and the probabilistic 
questions, like, e.g., what is the probability to obtain a 
certain knot type under given conditions (e.g., upon loop 
closure). Much less is known about the more physical 
aspects, which are how knots influence the properties of 
polymers. The simplest question to ask about physical 
properties is what the average spatial size is of a poly- 
mer loop whose knot type is quenched. To this end, J. 
des Cloizeaux .lj conjectured as early as 1981 that the 
size of a trivially knotted loop (i.e., an unknot) scales 
with the number of segments, N, in the same way as 
in the case of a self-avoiding walk, which is N v , where 
v = vsaw ~ 0.589 ~ 3/5. We should emphasize that 
the polymer in question is not phantom in the sense 
that segments cannot cross each other, but it is assumed 
to have a negligible excluded volume (or thickness) 0- 
Thus, according to des Cloizeaux's conjecture, exclusion 
of all knots acts effectively as volume exclusion. More 
systematic arguments, albeit still only scaling level, to 
support this conjecture were presented more recently in 
the work sj, yielding the following prediction for the 
(mean square) average gyration radius of a trivially knot- 
ted loop: 



(-Ro)triv 



(£ 2 /12) N if N < N 
A (£ 2 /l2) N 2v if N » N 



(1) 



Here, £ is the segment length, and the parameter Nq 
is sometimes called the characteristic length of random 
knotting; it appears in the probability of observing a triv- 
ially knotted conformation (an unknot) in a fluctuating 
phantom (i.e., freely crossing itself) loop: 



wtrw = w exp (-N/N ) 



(2) 



When N is smaller than Nq, a phantom loop has few 
conformations with non-trivial knots. Therefore, the set 
of allowed conformations for an unknotted non-phantom 
loop is not significantly different from that of a phan- 
tom loop. This is why at N < Nq the Gaussian scaling 
of gyration radius is expected. For this case, the £ 2 /12 
prefactor results from the facts that, (i), the mean square 
gyration radius for the linear chain is 1/6 of its mean 



square end-to-end distance, £ 2 N, and, (ii), (R 2 ) for the 
loop is half that for the linear chain Q. Prefactor A 
for the N > No regime in formula |QJ must provide for 
smooth cross-over between regimes at N ~ No, which 
means 



A{£ 2 jY£)N 2v k, (£ 2 / 12) N {N/N y 



A 



(3) 

Given the fundamental character of the problem, and 
given that theoretical arguments remain far short from 
mathematically rigorous, it is vital to look at the sim- 
ulation data. The situation on this front is at present 
contradictory. The difficulty is that the N" scaling is 
only expected at JV 3> Nq, while N , according to all 
simulations 0, |(| , although somewhat model-dependent 
(e.g., segments of fixed length vs. segments of Gaus- 
sian distributed length), is as large as around No ~ 300 
for some models. The work Q claimed a few data points 
consistent with the prediction , but its method of loop 
generation was later criticized y|. In the work the 
authors came to the contradictory conclusion that the 
N" scaling is observed upon fitting the R g dependence 
on N over the entire interval of N from well below to 
well above iVo, while in the N > No range the Gaussian 
behavior N 1 / 2 is recovered. This conclusion appears to 
suggest that the loops with N < No, which experience 
virtually no topological constraints, swell most strongly 
due to these constraints, which does not seem possible. 
Also, this result is not supp orted by the other earlier work 
from the same group, |10|. where authors mostly looked 
at the role of excluded volume, but also formulated the 
conclusion that " when N is large enough . . . the value 
of the exponent u/c" (for the given knot JC) "should be 
consistent with that of self-avoiding walks" . Finally, in 
the recent work |8( authors examined polymers of up to 
N — 600 segments, and claimed to observe the N v scal- 
ing. 

In all of the works 0, 0, , in order to extract the 
value of scaling exponent from the data, which (particu- 
larly in |8j) is almost entirely restricted to the cross-over 
range, authors fitted the data using the formula 



{R 2 ) = AN 2v 1 + B(N /N) A + 



(4) 



with adjustable parameters A, B, and v, usually assum- 
ing for simplicity A = 0.5. This approach, suggested in 
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is motivated by the analogy with the renormaliza- 
tion group treatment of the excluded volume problem. 
Unfortunately, this analogy itself hinges on the idea that 
the power v in formula is the same as that for self- 
avoiding walks, which is precisely the idea to be tested. 
Furthermore, formula ijljl. even if valid, is not the inter- 
polation working across the cross-over range from trivial 
to non-trivial scaling. Indeed, this formula does not yield 
Gaussian scaling N 1 / 2 in any range of N. 

In this paper, we systematically test the prediction QJ 
for the length up to N = 3000; this length is determined 
by our current computational capabilities, but it is also 
about the threshold above which excluded volume effects 
get significant for DNA 2]. Consistent with the pre- 
dicted value, we find v w 0.58 ± 0.02. Furthermore, we 
were able to examine the probability distribution of the 
gyration radius, and find, for instance, that trivial knots 
are noticeably less compressible than the average of all 
loops. 

The plan of our simulation is as follows. First, we 
generated loops of the length N divisible by 3 using the 
following method. To produce one loop, we generated 
7V/3 randomly oriented equilateral triangles of perime- 
ter 3£. We consider each triangle a triplet of head-to-tail 
connected vectors. Collecting all N vectors from N/3 tri- 
angles, we re-shuffled them, and connected them all to- 
gether, again in the head-to-tail manner, thus obtaining 
the desirable closed loop 0. For each loop, we compute 
the gyration radius 

1 N N 

^ = ^EE4> r« = l*Wjl . (5) 
i=i j=i 

Second, once a loop is generated, we determined its topol- 
ogy by computing the topological invariants. For the 
loops with N < 300, we used Alexander invariant A(— 1) 
and Vassiliev degree 2 and degree 3 invariants «2 and i>3 
[l4| . The loop was identified as a trivial knot when it 
yielded |A(— 1) = 1, w = 0, and v 3 — 0. For longer 
loops with N > 300, for reasons of computational impo- 
tence, we only used A(— 1) and «2 invariants, assigning 
trivial knot status to the loops with |A(— 1) = 1, and 
V2 = 0. The details of our computational implementa- 
tion of these invariants are described elsewhere 0] . Of 
course, because of the incomplete nature of topological 
invariants, our trivial knot assignment is only an approx- 
imation, and surely was sometimes in error. 

At every N, we continued generating loops until col- 
lecting the desirable number of presumably trivial knots, 
as specified in Table H] Collecting this amount of statis- 
tics required more than 10 5 CPU hours, roughly 11 CPU 
years. This extraordinarily long execution is the painful 
result of the exponentially rare nature of trivially knotted 
loops 0. 

The first result of our simulations, presented in figure 
n is the fraction of trivial knots among all loops, Wtnv, 
as it depends on TV. Overall, our data agree well with 
exponential formula @ and the data of earlier simula- 



TABLE I: Minimum number of loops used for statistics of 
each point 



N 


loops 


trivial knots 


CPU 




generated 


produced 


hours 


15 to 480 


10 6 


10 6 


0.02 to 288 


510 to 990 


7 x 10 s 


10 5 


41 to 1.7 x 10 3 


1020 to 1701 


7 x 10 5 


10 4 


230 to 8.2 x 10 3 


1800 to 2301 


10 6 


10 3 


3300 to 2 x 10 4 


2400 to 2502 


10 6 


10 2 


(3.6 to 7.7) x 10 3 


3000 


10 6 


9 


3.5 x 10 3 



fraction trivially knotted 




500 1000 1500 2000 2500 N 



FIG. 1: The fraction of loops generated with trivially knotted 
topology followed the well known exponential form as a func- 
tion of loop length N. Deviation from the fit line at large N is 
due to the incompleteness of topological invariants employed 
and reflects the contamination of the supposedly trivial pool 
with some non-trivial knots. The inset shows the same data 
for the interval of N up to 300, where third degree Vassiliev 
invariant was used along with second degree Vassiliev and 
Alexander invariants. 



tions |5j, |£j • However, deviation from the exponential is 
apparent in the region N > 1500. Later in this paper, 
we shall look more closely into this deviation. We be- 
lieve that it results from use of insufficiently powerful 
invariants in assigning trivial status to a loop, and re- 
flects the contamination of the supposedly trivial pool 
with some non-trivial knots. Accordingly, to extract the 
parameters Nq and wq, see (0, we used only data in the 
range 50 < N < 300, where the occurrence of mistak- 
enly identified knots is lower, and where we could rely on 
the third Vassiliev invariant in addition to the other two. 
This yields the best fit parameters Nq — 241 ± 0.6 and 
wo = 1.07 ±0.01. Our result for the characteristic length 
of random knotting is somewhat smaller than reported 
in the previous works 0, 0] , which we interpret as due to 
the subtle difference in the models examined |lfij |. 

We now approach the central issue of this paper, which 
is our data on the gyration radius of loops, as plotted 
in Figure El To begin with, as a consistency check, at 
each N we look at the (Rg) a \i averaged over all gener- 
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radius for all loops (which is equal to N£ 2 /12) reads 



all 



N£ 2 



= w(Rl) tTiv + (1 - w)(R 2 q ) non - triv . (6) 



(-Rg)triv - (^g)tr 



non — triv 



FIG. 2: Gyration radius averages over trivially knotted loops, 
and, as a control, over all loops. The trivial knot average 
exhibits power law behavior at large N similar to that ex- 
perienced by polymers which have excluded volume. This 
topology driven swelling is seen to develop beyond the crit- 
ical length No = 241. Independently collected data of 
is shown by stars (*) and agrees with our results. Upper 
Inset: Swelling parameter, (Rg) all/ (N£ 2 / 12), averaged over 
all loops irrespective of their topology, shows no dependence 
on N, lending credence to our loop statistics. By contrast, 
swelling parameter {R 2 g ) t rw / \N£ 2 /12) demonstrates that triv- coupled equations to find 
ial knots on average swell significantly above the Gaussian 
average size N£ 2 /12. Lower Inset: Trivially knotted gyra- (R 2 ) 
tion radius average normalized by power law fit to data. Solid '"' 
lines demarcate ±5% deviation of data from fit. The small 
difference between corrected (see formula JSJ), and raw trivial 
knot gyration radius suggests that errors in knot assignment 
do not significantly affect the scaling power v. 



Now let f5 be the probability of a non-trivial knot mis- 
takenly assigned as trivial. In Figure ^ 8 is visible as 
the vertical distance the data points rise above the fit 
line. The fraction of loops to which we assign, correctly 
or mistakenly, the trivial status is w + 5. The conditional 
probabilities of the loop to be a trivial or non-trivial 
knot provided it is assigned trivial status by our imper- 
fect topological invariants are w/(w + S) and S/(w + 5), 
respectively. Accordingly, the gyration radius averaged 
over thus contaminated trivial pool, (R 2 )' tliv , can be 
described as the weighted average of loops which cither 
posses or lack trivial knot topology, thus, 



(R 



non— triv " 



(7) 



Implicit here is an assumption that mistakenly iden- 
tified knots have the same average gyration radius 



as all other non-trivial knots. Accepting it, 



we observe that in the equations (JSJ and J7J) we know ev- 
erything except (i?^)triv and (-Rg)non-triv We solve these 



triv 



gl triv 



'g/triv 



j [(^g)triv ~ (-flg)all] 

w(l — w — S) 



(R 



1 



-) - (RlU (- 

w J y \w 



(8) 



ated loops, irrespective of their topology. As the upper 
inset of Figure [2] indicates, the swelling parameter, de- 
fined as (RT) a n / '(N I 2 1 '12) , is practically independent of 
N. Since N£ 2 /12, as we pointed out before, is the mean 
square gyration radius for Gaussian loops, this result con- 
firms the statistically representative character of our loop 
ensemble. 

Similar swelling parameter for trivial knots, 
(-Rg)triv / {N£ 2 /12) is shown in the same upper in- 
set of Figure [21 A few points independently collected by 
A.Vologodskii Qjj using only the Alexander invariant 
are shown as stars, they agree with our data. The data 
demonstrate clearly that loops with trivial knot topology 
are on average much more extended than other loops. 

To move beyond this qualitative conclusion to the 
quantitative characterization of topology-driven swelling, 
we found it necessary to look closer at the errors caused 
by the contamination of the trivial knots pool due to 
mistaken assignment of some non-trivial knots as trivial 
due to insufficiently powerful topological invariants. We 
used the following procedure to correct for this problem 
of trivial ensemble contamination. 

Let w = Wtiiv = wo exp(— N/No), the true probability 
of finding a trivial knot. Then, the averaged gyration 



where the later simplification makes use of the observa- 
tion that S <C w everywhere, and that w« 1 when the 
correction in question is noticeable (say, at N > 1000 
or so). Thus, we obtain that (R 2 )trW is somewhat 
larger than directly measured quantity (R 2 )' tliv (because 
(Rg)' tliv > (Rg) a ii) by the amount proportional to S/w. 

Thus corrected data for (-Rg)triv are presented in the 
main part of Figure [21 as diamonds (o) . The data fit 
well to the simple power law (ii|)triv = A (£ 2 /12) N 2l/ at 
N larger than about N — 500, with best fit parameters 
A w 0.44 ± 0.03 and v w 0.58 ± 0.02. This result is 
fully consistent with theoretical prediction 3] in several 
respects. First and foremost is the very fact of power 
law dependence of (-Rg)triv on N. Second, the value of 
exponent v matches well that of the self- avoiding walks, 
thus confirming the des Cloizeaux conjecture [l|. Third, 
the range of N where the power law is observed supports 
the idea that it should start at N > Nq, as in formula QJ. 
Fourth, the value of prefactor A agrees with prediction 
© which is A w 0.42, thus providing for smooth cross- 
over at N close to No, as expected. 

The fit quality is addressed in the lower inset of the 
Figure [21 where data/ fit is plotted against N. Over- 
all, data remain within ±5% of the fit. Importantly, the 
difference between corrected (see formula © ) and un- 
corrected data is within the 5% error corridor, thus sug- 
gesting that the fit result is reliable and is not affected 
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dramatically by the inevitable errors of knot identifica- 
tion. 

At the same time, we should point out that within the 
5% corridor, our data exhibit small but systematic bend 
upwards. Formally, this leads to the observation that 
power law fitting of only part of our data, starting from 
a larger N, say N > 1000 or N > 1500 etc, yields in- 
creasing v, up to the physically absurd values of 0.9 or 
so at very large N. Of course, these unphysical results 
come from narrowing windows of data where the statis- 
tics are increasingly poor. Nevertheless, currently we do 
not know if the upward bend of the data curve in Figure 
12 is entirely due to the measurement errors, or if it hints 
to something more serious. In particular, this bend pre- 
vents us from meaningfully fitting the data with formula 
Q. Further work is needed to understand whether data 
improvement, formula modification, or both is required. 

Within our current capabilities, we can use our data 
to address quite a few more interesting issues. One pos- 
sibility is to look at the mean square gyration radius of 
non-trivial knots. Such data are presented in Figure 
Apart from being pulled to much larger values of N, our 
data in this respect is quite similar to that presented ear- 
lier by the Swiss group 8J. For every non-trivial knot, 
the mean squared gyration radius remains smaller than 
the topology blind average over all loops, N£ 2 /12, and 
becomes larger at a certain value of N characteristic 
for each knot. On theoretical grounds, it was hypoth- 
esized [llj that the leading term in N — > oo asymptotics 
(R 2 g ) ~ A (£ 2 /12) N 2v should be valid for every particu- 
lar knot type, with both scaling power v and "amplitude" 
A independent of the knot type. Indeed, this conclu- 
sion seems inevitable considering the fact that any given 
knot at sufficiently large N is dominated by the stretches 
which effectively look like parts of a trivial knot (see also 
Looking at the data, Figure |3| with this theo- 
retical concept in mind, we see that the sizes of all knots 
considered do approach each other with increasing N. 
However, this happens quite slowly even when N is as 
large as, say, N — 1500 ~ 6N . 

Our data allow us to make one more step and to look 
not only at the averaged value of R 2 for trivial and some 
non-trivial knots, but also at the entire probability dis- 
tributions. We were able to generate and analyze his- 
tograms of quality (i.e. looking smooth when plotted) 
for loops of size N < 600, where contamination of the 
trivial pool and the corresponding correction © are to- 
tally insignificant. Predictably, the probability distribu- 
tions are different for different topological classes, such 
as all loops versus loops of a certain knot type K. Also 
predictably, the probability distributions of R 2 spread 
out as N increases. The latter observation suggests the 
idea of looking at the probability distributions of the re- 
scalcd variable p = R 2 /(R 2 ), where the normalization 
factor (R 2 ) is taken separately for each N and for each 
topological entity. 

Our main findings are summarized in Figure 0] where 
we present probability distributions P(p) for the trivial 



<R a 2 > K /(NI 2 /12) 





1000 2000 



50 100 200 500 1000 2000 N 

FIG. 3: Log-log plot of the mean square gyration radius, 
(Rg) K , of knot type «, normalized by the topology blind aver- 
age over all loops for several particular knot types. The inset, 
which shows the ratio of a particular knot gyration radius to 
the trivial knot gyration radius, (Rg)o, demonstrates that all 
knots remain smaller than, but approach the size of trivial 
knots. 



knots 0i (o), trefoils 3i (A), and 4i knots (□). The data 
shown are for N = 90, where high quality statistics were 
most easily attainable: each histogram is the result of 
more than 20 million loops. 

In the same figure 21 we plot also for comparison the 
analytically computed probability distributions for linear 
chains and for all loops. For linear chains, the necessary 
distribution Pii nea r(p) was found by Fixman a long time 
ago [19]. He showed that the corresponding character- 
istic function (Fourier transform of probability density) 

is equal to K\i neax (s) = (sin z/z) 3 ^ 2 , where z 2 — Ais, 
and where s is conjugate to p (i.e., Fourier transform in- 
volves e lsp ). We were able to derive similar expression 
for the probability distribution over all loops, irrespec- 
tive of topology. In this case, the characteristic function 
reads if all loops (s) = (2 sin(,2/2)/z) _3 , where z 2 = 8ts, 
with the same definition of s. Numerical inversion of 
Fourier transforms yield the curves presented in Figure 
0] To avoid overloading the figure, we do not show the 
corresponding data points obtained for linear chains and 
for all loops, but they all sit essentially on top of the the- 
oretical curves (which is comforting, as it confirms once 
again the ergodicity of our loop generation routine) . 

Comparing the shapes of probability distributions for 
all loops and those with identified topology, we notice 
that the latter distributions are somewhat more narrow. 
Although the effect looks small for the eye, it is certainly 
there, and it is well above the error bars of our measure- 
ments. This means simple knots are less likely to swell 
much above their average size than other knots, and they 
are also less likely to shrink far below their average, again 
compared to other knots. 

The latter point is of particular interest given its re- 
lation to all problems involving collapsed polymers, such 
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FIG. 4: Probability density plots for chains (line, [l9jp. all 
loops (another line), and loops with certain knots (Oi - o, 3i - 
A, 4i - □). Distributions are presented in terms of the scaling 
variable p = R^/(R? g ). Lower Inset: Semi- log probability 
density plot (or linear entropy plot) at large p. Upper Inset: 
Semi-log probability density plot (or linear entropy plot) at 
small p against 1/p. 



as proteins. Closer look at the small R g region of the 
probability distribution is presented in the upper inset of 
Figure 0] There, the probability distributions are plot- 
ted in the semi-log scale against 1/p. This can be also 
understood as the plot of " confinement" entropy, which 
corresponds to the squeezing the polymer to within cer- 
tain (small) radius. The reason why we plot the data 
against 1/p is because both Pn ne ar(p) and P a ii loops at 
small p have asymptotics ~ exp (— const/p), which corre- 
sponds to confinement entropy ~ 1/p, and which can be 
established by a simple scaling argument, as described, 
e.g., in ,4'] (page 42). This 1/p behavior is seen clearly in 
the upper inset in Figure 0] Furthermore, we see indeed 
that compressing any specific knot, trivial or otherwise, 
is significantly more difficult than compressing a phan- 
tom loop. Analytical expression of entropy for knots is 
not known, only the R~ 3 ~ p~ 3 / 2 scaling at small p was 
conjectured in the work [20|. Although our data is qual- 
itatively consistent with this prediction in terms of the 
direction of the trend, more data is needed for quantita- 
tive conclusion. 

To conclude, we want to speculate on some broader im- 
plications of our findings. To this end, it seems obvious 
that for any given knot IC, the final asymptotics of (ic)jc 
at N — > oo are governed by the same exponent v as for 
the trivial knots, that is, according to our results above, 
Vfz = vtiiv — ^saw ~ 0.589 w 3/5. To explain this point, 
and also to understand at which N this asymptotics takes 
over, it is useful to compare the polymers with quenched 
and annealed topology, the latter being simply phantom. 
At every N, the polymer with annealed topology samples 
a certain ensemble of knots; as N grows longer, the set 
of knots which are sampled becomes more diverse. In a 
loose sense, we can imagine a certain average for this set 
of knots, something like average number of knots, or aver- 



age knot complexity. For instance, we can average diame- 
ter of maximally inflated tube |2l| or minimal rope length 
|22| . Let us denote IC* (N) as the average, or typical knot 
for the given length N. Clearly, as N increases the typi- 
cal knot gets more complex, its rope length increases, and 
its inflated diameter decreases. Now let us go back and 
consider the real polymer of the given length N with real 
quenched knot topology, IC. We should recognize the im- 
portant difference between the cases when given knot /C 
is more complex (with longer rope length or smaller tube 
diameter) or simpler (with shorter rope length or larger 
tube diameter) than IC*(N). In the former case, our real 
polymer can be called overknotted, because it contains 
larger amount of knots than it would do spontaneously 
if allowed. In the latter case, the polymer can be called 
underknotted, because it has fewer knots than typical for 
its length. Overknotted polymers should be more com- 
pact compared to the annealed, or phantom loop; in other 
words, for them we expect {R 2 9 )k < N£ 2 /12. By contrast, 
underknotted polymers should be more swollen compared 
to their phantom counterparts, (R 2 )k > N£ 2 /12. In the 
light of this consideration, we can now understand what 
happens if we have a given non-trivial knot /C and we 
increase N. At the beginning, N is small and we are in 
the overknotted regime. Eventually with growing N we 
expect to cross over into the underknotted regime, and 
it is in this regime that we expect the size of the knot 
to scale as N", because every underknotted loop consists 
basically of very long pieces which are not entangled with 
each other and are not knotted themselves. The number 
of such pieces depends on the knot /C, but does not de- 
pend on N, such that their length scales as N and their 
size, therefore, must scale as N v . Much work is needed 
to make these considerations more quantitative and less 
speculative, a start in this direction would be to define 
K.*{N) in a more rigorous fashion. Among other things 
relevant here, the issue of knots localization [23l l24l |25| 
must be quantified and incorporated. 

To summarize, we have presented simulation data on 
the sizes of loops with the topologies of trivial or non- 
trivial knots, for the lengths of up to 3000 segments. We 
found that topological constraints have marginal effect on 
the loop size as long as the loop is shorter than the char- 
acteristic length of random knotting, which is about 250. 
At larger N, our results for trivial knots are consistent 
with crossing over into the scaling regime R g ~ N" analo- 
gous to self-avoiding walks statistics, for which v « 0.59. 
Our findings are also consistent with the idea that the 
size of any particular non-trivial knot becomes asymp- 
totically equal to that of the trivial knot at very large N , 
although our data suggest the slow decaying approach to 
this asymptotic regime. Finally, looking at the proba- 
bility distributions of the (properly re-scaled) loop sizes, 
we found that topologically complex loops are less likely 
to adopt either strongly collapsed or strongly expanded 
configurations. 

We wish to thank A. Vologodskii for sharing with us 
his unpublished data ^?|]- We also acknowledge fruitful 
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discussion with T. Deguchi. 
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